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We point out that in extended supersymmetric models such as supersymmetric left-right models, 
it is possible to have leptoquarks that explain the HERA high Q 2 anomaly without giving up 
R-parity conservation. The leptoquarks belong to vectorlike (2, 2, ±|, 3or3*) representations of 
SU(2)l x SU(2)r x U(1)b-l X SU(3) C (denoted by G2213)- Considerations of D-terms imply that 
the only phenomenologically viable model is the one where the leptoquark couples to positrons 
and the up quark. Unlike the R-parity violating scenario, the leptoquarks are accompanied by 
new superpartners,the leptoquarkino which leads to many intersting signatures in other collider 
experiments. At Tevatron, pair productions of the leptoquarkino will give rise to dilepton signals 
very distinct from the top productions. These models can lead to unification of gauge coupling 
constants at a scale of around 10 10 GeV implying that grand unification group is not of the usual 
SU(5) or 50(10) types but rather an automatically R-parity conserving 517(5) x 5(7(5) GUT model 
recently proposed by one of the authors (R. N. M.) which leads to a stable proton. 
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If the high Q 2 anomaly observed recently in the e + p 
scattering by the HI [jlj and the ZEUS collaborations 
is confirmed by future data, it will be an extremely inter- 
esting signal of new physics beyond the standard model. 
A very plausible and widely discussed interpretation of 
this anomaly appears to be in terms of new scalar par- 
ticles capable of coupling to e + u or e + d of scalar lepto- 
quarks H with mass around 200 GeV. Alternative inter- 
pretations based on contact interactions Q or a second 
Z' H have been proposed; but attempts to construct 
models that lead to the desired properties seem to run 
into theoretical problems. 

The leptoquark must be a spin zero color triplet par- 
ticle with electric charge 5/3 or 2/3. The latter electric 
charge assignment makes it possible to give a plausible 
interpretation of the leptoquark as being the superpart- 
ner of the up-like quark |j) of the minimal supersym- 
metric standard model (MSSM) provided one includes 
the R-parity violating couplings of the type \'QLd c to 
the MSSM (where Q, L denote the quark and the lepton 
SU(2) L doublets and u c ,d c ,e c denote the SU(2) L sin- 
glets [Q ) . There are however very stringent upper limits 
on several R-parity violating couplings: for instance, if in 
addition to the A' term described above, one adds the al- 
lowed \"u c d c d c terms to the superpotential, then it leads 
to catastrophic proton decay unless A' A" < 10~ 24 [||. 
There are also stringent limits on X' nl < 10 -4 |J from 
neutrinoless double beta decay, which forces the lepto- 
quark to be c or t rather than the obvious choice u. Fur- 
thermore, within such a framework, the lightest super- 
symmetric particle (LSP) is no more stable and there- 
fore, there is no cold dark matter (CDM) candidate in 
such theories. In view of the fact providing a natural 



CDM candidate was long considered an attractive feature 
of the supersymmetric models, it may be worthwhile to 
consider extensions of the MSSM which incorporate the 
leptoquark without giving up R-parity conservation (and 
hence the idea of LSP as a natural CDM candidate) . It 
is the goal of this letter to report on the results of such 
a study. 

The class of supersymmetric theories where R-parity 
conservation is automatic can provide an absolutely sta- 
ble LSP that can act as the CDM candidate. We will 
therefore use them as a typical framework for studying 
the consequences of leptoquarks incorporated into such 
theories while at the same time maintaining automatic 
R-parity conservation. Minimal versions of such theo- 
ries are based on the gauge group SU(2)r, x U(l)i 3R x 
[7(1) b-l or SU{2) L x SU{2) R x U{1) B -l © with usual 
assignments for the quarks and leptons. We will con- 
sider the latter left-right symmetric gauge group. The 
quark doublets Q = (u, d) and Q c = (u c , d c ) transform as 
doublets under the SU(2)l and SU{2)r groups respec- 
tively and similarly the lepton doublets L = (y, e) and 
L c = {v c ,e c ) respectively. The SU{2)r gauge symmetry 
breaking is achieved via the B — L non-singlet isotriplets 
A c ee (1, 3, -2) and A c = (1, 3 + 2) and their lefthanded 
counterparts added to maintain left-right symmetry (the 
numbers in the parenthesis denote the SU(2)l, SU(2)r 
and U(1)b—l quantum numbers). The standard model 
symmetry is broken by the bi-doublet fields 0(2, 2, 0). 

We augment this model by including the leptoquark 
fields. Demanding that the lepto-quarks couple to e + d 
or e + u leads to the conclusion that they must be- 
long to the multiplet (2,2,4/3,3) (denoted £ql<0 and 
(2, 2, -4/3, 3*) (denoted by T,q-l)- Each of these multi- 
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plets have four scalar fields and four fermion fields which 
will be denoted in what follows by the obvious subscript 
corresponding to their couplings. We denote the four 
scalar leptoquarks as £ ue c , E u!/ c , £d e c and £di/<= and their 
fcrmionic partners (to be denoted by a tilde on the cor- 
responding scalar field). Writing down the most general 
superpotential, one can easily convince oneself that the 
resulting theory maintains the property of automatic Im- 
parity conservation. Thus the lightest neutralino LSP 
will be stable in this model and can serve as the cold 
dark matter. 

Before discussing the application of this model to ex- 
plain the HERA anomaly, let us first discuss the mass 
spectrum of the model. The superpotential for this model 
will have a direct mass term of the form Mo££ which will 
imply that the fermionic fields in the leptoquark multi- 
plet will have a masses Mo prior to symmetry breaking. 
There may be other contributions to these masses from 
radiative corrections which will split the degeneracy im- 
plied by the above mass term. 

As far as the scalar leptoquark states are concerned, 
their masses will receive several contributions: first a 
direct common contribution from the Mo term given 
above. After symmetry breaking, the D-terms of the 
various gauge groups will contribute. There is also soft 
SUSY breaking contribution along with the radiative cor- 
rection. Assuming that the SU(2)r symmetry is bro- 
ken by the vev's < A c >= V\, < A c >= v-i and 
the SU(2)l symmetry is broken by the two <p vev's as 
diag < <p u > = (0, vsin/3) and diag < <fid >— (vcos0,O), 
we can write the masses for the various scalar leptoquarks 
as follows: 

Ml = M 2 + {I a R gl R - ?-±gl_ L ){vl - vl) (1) 
+ 92lV IsL cos2(3 + + Radiative correction 



A^ is the soft SUSY breaking contribution. The values 
of the If, and B — L are given in Table I: 
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TABLE I. The I^r and Bl quantum numbers of the vari- 
ous leptoquark states 



The radiative correction can be positive or negative. 
From table 1 and Eq.(l), we see easily that if V\ < V2, 
then the lightest leptoquark state is the first entry in 
the table I which corresponds to the conjugate of T> ue c 
since g\ R > 2g 2 B _ L in the left-right models. Furthermore, 
interestingly enough for this choice of vev's, assuming the 
combined A^ + Radiative correction to be smaller than 
the the D-term, the leptoquarkino states are heavier than 
the lightest leptoquark state. We assume their masses to 
be in the range of 300 to 400 GeV. In a subsequent section 
we will obtain a lower bound on the leptoquarkino mass 
from the present collider data. 

Turning to the couplings of the leptoquarks £ and £ 
to quarks and leptons, it is given by the superpotential: 

Wi q =\ ij {m c i Lj+VQiL<i) (2) 

(where i, j are generation indices). Let us assume for sim- 
plicity that Xij are diagonal in the mass eigenstate basis 
for the quarks and leptons. Explanation of the HERA 
high Q 2 anomaly seems to require An ~ 0.05 which we 
will assume from now on and mass of the scalar field 
in £ (assumed to be £ u c e from the above mass argu- 
ments) around 200 GeV. As far as the other couplings 
go, they will be strongly constrained by the present ex- 
perimental upper limits on the low energy processes such 
as /i — > e7, r — > e7 etc. which will arise at the one loop 
level from the exchange of £ and £. There are also tree 
level diagrams which can lead to rare processes such as 
K — > 7re~/i + . These processes imply an upper limit of 
A22 < 2 x 10~ 3 . The most stringent upper limit on A22 
comes from the upper limit on the process — > /i + e~ 
and yields A22 < 2 x 10~ 4,5 . Turning to A33, the most 
stringent limits arise from the present upper limit on the 
branching ratio for the process r — > e"f which the 1996 
Particle data tables give as < 1.1 x 10~ 4 [jll]]. This im- 
plies a weak upper limit on A33 < 0.2. Thus the third 
generation leptoquark coupling A33 could in principle be 
comparable to An. If the efficiency for the detection of 
r leptons at HERA were comparable to the detection ef- 
ficiency for electrons, that could also sever ly limit A33. 
Finally we also note that such a value for An is also con- 
sistent with the present data on the parity violation in 
atomic physics. Thus it appears that our leptoquark cou- 
plings are consistent with all known low energy data. It 
is clear from the above discussion that the leptoquarks 
do not respect [i — e universality. 

The existence and properties of the leptoquarkino pro- 
vides a new and unique signature for our model as 
compared to all other proposals to explain the HERA 
anomaly. To see this note that in hadron colliders, we 
can produce pairs of leptoquarkinos at the same rate as 
the ti pair due to identical color content. Moreover, due 
to R-parity conservation the leptoquarkino decay leads 
to a missing energy signal as follows: £ — > e + u with 
ii — > 11 or £ — > u -\- with e+ — > e + + x ■ In both 
the cases we have e + , u plus missing energy in the final 
state. If the A33 coupling is comparable to An as is al- 
lowed, then the branching ratio for leptoquarkino decay 
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to electrons will be 50%. This signal is similar to the 
top signal at the Tevatron with one crucial difference. 
The dilepton branching ratio from the leptoquark pairs 
is 100% compared to only 10% from the top pairs. There 
will be no events. The branching ratio to e + e~ 

will be 25% compared to only 1% from the top. The 
present observations should therefore lead to lower limits 
on the mass of the lightest leptoquarkino pairs. Thus an 
excess of dilepton pairs over that expected from the top 
productions, or an excess in e + e~ channel over fi + fj,~ will 
be a clear signal of leptoquark productions at Tevatron. 
In addition, leptoquarks will also give rise to harder lep- 
tons and larger missing energy events. In Fig. 1, we plot 
the cross section for the leptoquarkino pair production. 
We see that for the present combined Tevatron data of 
about 200 pb _1 , there would be 6 events of type e + e~jj 
plus missing energy for Mf, c = 250 GeV, assuming a 
detection efficiency of 0.2 for the e + e~ mode. At present 
CDF collaboration has only one such event jlj] , and the 
DO collaboration has also 1 event [|lj|. Using these two 
events, and the the cross sections given in Fig. 1, we ob- 
tain a lower bound on the leptoquarkino mass of about 
290 GeV ( assuming the detection efficiency to remain 
0.2 for the higher masses). Any evidence for an excess of 
dielectron events would be interpretable interms of this 
new leptoquarkino. We urge the CDF and DO colubo- 
ration to look for such excess, and also to look for any 
dilepton event in which the lepton Pt or missing energy 
do not fit the top productions. The recent CDF analysis 
shows that the leptoquark mass bound at 95% CL 
is 210 GeV assuming the branching ratio is 1. In this 
model however the branching ratio is less than 1, since 
the leptoquark can also decay into the top quark and the 
r lepton. 

So far we have discussed the signals of leptoquarkinos 
in the usual SUSY theories. If however the supersymme- 
try breaking is communicated to the obseravable sector 
by the SM gauge ineteraction the signal is differ- 
ent. Let us first consider the case when neutralino is the 
NLSP in these models. The leptoquarkino decays into a 
spositron and u quark (squarks masses are large in these 
models). The spositron decays (100%) into a positron 
and a neutralino and the neutralino then decays (100%) 
into a hard photon and a gravitino. Consequently the 
final state in the leptoquarkino pair production process 
has e + e~77jj plus missing energy without any SM back- 
ground. If lighter stau is the NLSP and neutralino is the 
NNLSP, spositron decays (100%) into a positron and a 
neutralino and the neutralino then decays (100%) into a 
r and a lighter stau. The lighter stau then decays into a 
tau and a gravitino. Consequently the final state in the 
leptoquarkino pair production process in the Tevatron 
has e + e~2T + 2T~jj plus missing energy. Out of these six 
leptons in the final state, one r pair (produced from the 
decay of stau) has much higher Pt than the other lep- 
tons. 

One could also look for the signals of leptoquarkino 



production in e + 7 colliders. The leptoquark or lep- 
toquarkino can be singly produced in such a collider. 
The scalar leptoquark will be produced in the associ- 
ation of an antiquark |lq| . As discussed before in the 
usual SUSY theories, the leptoquark further decays into 
a quark and a positron with the final state consisting of 
positron +jets. The leptoquarkino however will be pro- 
duced along with a squark. This leptoquarkino will then 
decay into a lepton and a squark. The final state has 
electron +jets+missingenergy. The missing energy part 
will then disentangle the leptoquarkino signal. In the 
gauge mediated SUSY breaking scenario, the final state 
has either a hard photon or t + t~ along with electron 
+jets+missing energy. 

In e + e~ collider the signal for this leptoquarkino would 
be jets -(-missing energy, since a squark anti-squark pair 
will be produced from a t channel exchange of lepto- 
quarkino. The exchange of a t channel leptoquark how- 
ever gives rise to jets without any missing energy. It is 
not possible for LEPII to see the sign al of leptoquarkino 
due to the large squark mass. A higher energy machine is 
needed for that purpose. A detailed study of the discov- 
ery prospects of the leptoquarkino in different colliders 
will be presented elsewhere. 

Another important implication of our model is that 
there does not seem to any operator which can give a 
non-negligible charged current signal at HERA. Thus ob- 
servation of charged current like events above background 
at HERA will be an evidence against this model. 

Let us now discuss the possible implications of the ex- 
istence of a low mass leptoquark supermultiplet for uni- 
fication. We will assume that the physics immediately 
beyond MSSM (i.e. in the TeV region) is described by 
a supersymmetric left-right symmetric model (SUSYLR) 
with the addition of the leptoquark multiplets described 
above. The evolution of the gauge couplings depend not 
only on the individual beta functions but also on the 
nature of the final unification group which determines 
the normalizations of the various gauge couplings at low 
energies. For instance, if we envision embedding the SU- 
SYLR group within a simple group such as 50(10) or 
SU(16) etc, the properly normalized weak hypercharge, 

Iy is given by the familiar formula ly = and 
the properly normalized B — L charge Ibl is given by 
Ibl = \[^{ B ~2 L )- This leads to the matching formula for 

the weak hypercharge coupling to be aj = fct^R + la^ 
at Mr. Evolving our model with this kind of unification 
leads to a unification scale of about 10 9 5 GeV assum- 
ing that at the Wr scale, we have two bidoublets, right 
handed triplets of type (1, 3, —2, 1) + (1, 3, +2, 1) without 
their left-handed partner in addition to the aforemen- 
tioned leptoquark multiplets. This is therefore unaccept- 
able since it will lead to catastrophic proton decay via the 
exchange of gauge bosons that violate baryon number in 
the 50(10) models. We therefore consider a different 
kind of embedding of SUSYLR into an SU(5) x SU(5) 
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model suggested in a recent paper p7[ . 



SU(5) x SU(b) embedding of SUSYLR with leptoquarks 



First let us briefly recapitulate the fermion assign- 
ments of the model: they belong to the (5 + 10, 1) + 
(1,5 + 10) multiplet JlT]] for every generation. One 
therefore needs a weak singlet vector like pair of D, U 
quarks and a heavy weak singlet vectorlike charged lep- 
ton E ± . The lefthanded (5,1) multiplet then consists 
of (D\, D%, -Dg, e~, v e ) whereas the righthanded multiplet 
(1, 5) is given by {D\, D2, D 3 , e c , v c ). The assignments to 
the 10 dimensional representations are easily obtained 
jl7| . For instance, the fermion assignment in (10, 1) rep- 
resentation is given by 
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and similarly for the (1, 10) multiplet which contains the 
other chirality states for the above fields arranged ex- 
actly the same way. The Higgs mechanism of the model 
is implemented by multiplets of type (5, 5) + (5, 5) and 
they lead to the vectorlike quarks and leptons acquiring 
mass at the unification scale. This leaves the low energy 
theory to be the usual SUSYLR model. The SUSYLR 
symmetry is broken down to MSSM by the Higgs mul- 
tiplets of type (1,15) + (15,1) which also leads to the 
see-saw mechanism for neutrino masses. Below the Wr 
scale the theory is MSSM with the important difference 
that now the R-parity violating couplings are automati- 
cally absent. In order to accomodate the leptoquarks in 
this model, we include the multiplets (5,10) + (10,5). 
We assume that below the GUT scale, only members of 
these multiplets that remain light are the leptoquarks 
discused above. 

Turning to unification with leptoquarks in this model, 
we first note that since U{\)b~l and SU(3) C now emerge 
from two different SU(5)'s and that there are more 
fermions in the fundamental representation of the GUT 
group than the 50(10) case, the normalization of the 
low energy couplings are totally different from the pre- 
vious case. For instance, now Iy — J i%-(y) an d Ibl = 
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\{IaL + I an) where a 
L generators. Because of this 



10 v 2 
denotes SU(S) 
unification profile is now very different; however with the 
same field content at Mr as above, we find a unification 
scale of 10 9 GeV in the one loop approximation. In Table 
1, we have showed some unification scenarios by varying 
the numbers of bidoublets. We have one (2,2,4/3,3) + 
(2, 2, -4/3.3*) and one (1, 3, -2, 1) + (1, 3, 2, 1) and have 



the Mz scale. From the table I, one can see that one can 
have one unified coupling at the GUT scale using seven 
bidoublets. It is important to point out that this low 
unification scale does not conflict with proton lifetime 
since proton is absolutely stable in this model |L7]]. The 
essential reason for this is that the vectorlike quarks and 
the light quarks do not mix with each other due to the 
choice of the gauge group and Higgs representations. 

We find it interesting that CDF collaboration, in their 
analysis of the dilepton channel, have reported |ll| 4 
events containing at least one r (3 of the events have 
at least one b quark). Their estimated background from 
non-tt productions is 2 events. Using their measured 
cross-sections of 7.5 pb (from DIL, SVX and SLT chan- 
nels), and their r channel acceptance of 0.12%, only one 
such event is expected from the tt productions in the 110 
pb _ 1 data. It is entirely possible that this event is due 
to statistical fluctuation, or due to charged Higgs pro- 
duction, however such excess events are expected from 
the leptoquark or leptoquarkino pair productions. One 
leptoquark decays to a r and the top quark, and the 
other to a e or a r. A leptoquark (or leptoquarkino) pair 
production will have a branching ratio of about 0.75 for 
the dilepton channel with a r in our model (compared to 
about 0.05 in tt productions). This will give an accep- 
tance (branching ratio x detection efficiency) of about 
1.8 %. Thus leptoquark (or leptoquarkino) pair produc- 
tion cross-section of about 0.5 pb will give rise to such 
an excess event. As can be seen from fig 1, this corre- 
sponds to produced leptoquarkino mass about 250 GeV 
(or a leptoquark of somewhat lower mass). 

In conclusion, we have suggested an alternative lep- 
toquark interpretation of the HERA high Q 2 anomaly 
within an extension of the MSSM in such a way that 
R-parity conservation is maintained automatically. The 
superpartner of the leptoquark, the leptoquarkino in this 
model leads to many interesting and testable predictions 
in present and future hadronic as well as in the future 
e-gamma colliders. At Tevatron, the pair productions of 
these leptoquarkinos give rise to very distinctive signal, 
namely opposite sign dilepton events of the type ee, er, 
rr accompanied by dijets plus missing energy, but no 
events of fifi type. If the leptoquarkino is not too much 
heavier than the "HERA leptoquark" , then these events 
will be observable in the present or upgraded Tevatron. 
It is possible that one excess event in the dilepton chan- 
nel containing a t, as reported by the CDF collaboration, 
could be due to the production of a leptoquark or lepto- 
quarkino. Although clearly more data is needed before 
definitive conclusion can be drawn. 
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TABLE II. The Grand unification scale, the intermediate 
scale and the strength of the unified couplings at the grand 
unified scale are shown for different numbers of bidoublets. 
ac A is the unified coupling for one SU(5) group and aa B is 
the unified coupling for the other 5(7(5) group 
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FIG. 1. Cross section for the production of leptoquarkino 
pair at Tevatron energy against the leptoquarkino mass. 
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